arXiv: 1506.06955vl [cond-mat.soft] 23 Jun 2015 


General contact mechanics theory for randomly rough surfaces with application to 

rubber friction 

M. Scaraggi 

Dll, Universitd del Salento, 73100, Italy and 
PGl, FZ Jiilich, 52425 Julich, Germany 

B.N.J. Persson 

PGl, FZ Julich, 52425 Julich, Germany and 
www.MultiscaleGonsulting. com 

We generalize the Persson contact mechanics and rubber friction theory to the case where both 
surfaces have surface roughness. The solids can be rigid, elastic or viscoelastic, and can be homoge¬ 
neous or layered. We calculate the contact area, the viscoelastic contribution to the friction force, 
and the average interfacial separation as a function of the sliding speed and the nominal contact 
pressure. We illustrate the theory with numerical results for a rubber block sliding on a road sur¬ 
face. We find that with increasing sliding speed, the influence of the roughness on the rubber block 
decreases, and for typical sliding speeds involved in tire dynamics it can be neglected. 


1 Introduction 

Practically all surfaces in Nature and technology have 
surface roughness over many decades in length scales, of¬ 
ten extending from the macroscopic size of the system 
down to atomic distances [l|. For this reason the inter¬ 
action between bodies in relative motion is a very com¬ 
plex topic. Nevertheless, because of its crucial impor¬ 
tance for many important tribological problems, such as 
the electric and thermal contact resistanceja-Q , rubber 
friction[l-Q, sealing Hi, frictional heating Q and tribo- 
electrification [l(i [lll| . a large effort has recently been de¬ 
voted to the contact mechanics between solids with rough 
surfaces d, [T2l - [l^ . 

Archardf^. and Greenwood and Williamson 
(GW)[2l[, have presented pioneering studies of the 
contact between solids with roughness. These theories 
approximate the surface asperities with spherical cups 
to which they apply the Hertz contact mechanism. The 
theory of Archard is for a rather idealized situation, 
so most contact mechanics studies have used the GW 
theory. However, the GW theory neglects the elastic 
coupling between the different contact regions (i.e. the 
lateral or horizontal coupling), as well as the vertical 
coupling between hierarchically distributed asperities 
(small bumps on top of bigger bumps on top of bigger 
bumps, and so on). For surfaces with roughness on many 
length scales this approximation leads to qualitative 
wrong results both for the contact area and the average 
interfacial separation as a function of the applied normal 

force, and other properties 

A very different approach was developed by Persson^ 
in 2001. This theory includes the elastic coupling and can 
be applied for arbitrary large normal forces, e.g., even as 
the contact approaches full contact, which is very impor¬ 
tant for some applications like the leak-rate of seals. The 
theory is very flexible and can easily be applied to com¬ 
plex situations, e.g., layered materials and elastoplastic 
or viscoelastic materials. The original theory was devel¬ 
oped for a viscoelastic solid with a perfectly flat surface 


in contact with a hard, randomly rough, substrate. In 
this paper we extend this theory to the case where both 
solids have surface roughness, and arbitrary viscoelastic 
properties, e.g., a viscoelastic solid in contact with an 
elastic solid. 

Consider the frictionless contact between elastic solids 
with nominally flat surfaces in normal (quasi-static) ap¬ 
proach. Assume that the solids have roughness described 
by the height profiles /ii(x) and /i 2 (x) [where x = {x,y) 
is a two-dimensional position vector in the surface xy- 
plane], and different elastic properties (Young modulus 
El and E 2 , and Poisson ratio ui and 1 / 2 , respectively). 
If the root-mean-square slope of the roughness on both 
surfaces is small one can show that this contact problem 
can be mapped on a simpler problem with one elastic 
half space with a perfectly flat surface, and another rigid 
solid with the combined roughness (see Fig. ??)[2^: 

/i(x) = /i 2 (x) — ft.i(x). 

The effective elastic module E and Poisson ratio iz of the 
elastic solid must satisfy 

1 — 1 — i^i 1 — t'l 

E ~ El ^ E 2 

We will refer to this transformation as surface-mapping. 
For viscoelastic solids both the Young’s modulus and the 
Poisson ratio depends on the frequency uj (and hence, 
during sliding (velocity v), on the length scale v/ui). In 
this case the surface-mapping is no longer valid, and in 
this paper we will show, within the Persson approach, 
how to treat the more general problem of viscoelastic 
contact mechanics when both solids have surface rough¬ 
ness. 

2 Qualitative discussion 

Let us consider a rubber block sliding on a rigid solid 
substrate. We assume first that the rubber block has 
a perfectly flat surface and the rigid substrate a ran- 
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Figure 1. Two solids with the surface profiles hi(x) and 
/i2(x, t)- The gap /i(x) = h2(x,t) — fii(x). If the upper solid 
moves at a constant velocity v = (vx,Vy) parallel to the bot¬ 
tom solid, we have h 2 (x, t) = /i 2 (x — vt). 


Figure 2. (a) A rubber block with a smooth surface sliding 
relative to a hard, randomly rough, substrate, (b) A rubber 
block with a rough surface sliding relative to a hard solid with 
a fiat surface. 


domly rough surface [see Fig. ??(a)]. In this case during 
sliding the asperities of the rigid solid will induce time- 
dependent deformations of the rubber, which results in 
the conversion of translation energy into heat via the 
internal damping in the rubber. Thus the viscoelastic 
deformations of the rubber will contribute to the kinetic 
friction coefficient. 

Let us now consider the opposite case where the rigid 
solid has a perfectly smooth surface and the rubber block 
a rough surface [see Fig. ??(b)]. In this case during 
sliding the deformations of the rubber does not change in 
time and no viscoelastic energy dissipation occurs. Thus 
in this case there is no contribution to the friction from 
viscoelastic deformations of the rubber. When roughness 
occurs on both surfaces the situation is more complex but 
already the argument given here indicates that the energy 
dissipation is dominated by the roughness in the harder 
solid. 

Another important difference between the two cases in 
Fig. ?? relates to the area of real contact. In Fig. ??(a) 
the contact area depends on the viscoelastic modulus of 
the rubber at finite frequencies. That is, during sliding 
at the speed v a surface roughness component with wave¬ 
length A will result in oscillating deformations of the rub¬ 
ber surface with a frequency a; ~ u/A. Since the elastic 
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Figure 3. If the contact area is small compared to the nominal 
contact area, introducing a periodic roughness on the rubber 
block result in a negligible change in the friction coefficient. 
Here we have assumed that adhesion does not manifest itself 
macroscopically as a pull-off force. 


modulus of rubber-like materials may increase by a factor 
of ~ 1000 when the perturbing frequency increases from 
the rubbery region to the glassy region, it is clear that a 
very strong dependency of the contact area is expected: 
An increase in the elastic modulus by a factor of ^ 1000 
can result in a decrease in the contact area by a factor 
of ^ 1000. However, when the surface roughness occurs 
on the rubber surface while the hard counter-surface is 
perfectly flat, there are no time-dependent viscoelastic 
deformations of the rubber and the contact area is deter¬ 
mined by the low-frequency viscoelastic modulus. 

Let us now compare the contact between a randomly 
rough hard solid (the substrate) and (a) a rubber block 
with a flat surface [see Fig. ??(a)], and (b) with a rub¬ 
ber block with a simple periodic “roughness” as indicated 
in Fig. ??(b). We assume that the longest wavelength 
roughness component of the substrate surface is smaller 
than the linear size of the (nominal) contact regions in 
??(b). A uniform stress or pressure cto is assumed to 
act on the upper surface of the rubber block. In case 
(a) a similar pressure is assumed to act at the interface 
between the block and the substrate}?^. According to 
the Persson contact mechanics theory, if frictional heat¬ 
ing can be neglected, the viscoelastic contribution to the 
friction for small applied load is proportional to the nor¬ 
mal load, i.e. the friction coefficient is independent of 
the load. Since the load is the same in both Fig. ??(a) 
and (b) it follows that the friction coefficient is the same 
too. This is only true as long as adhesion does not mani¬ 
fest itself macroscopically as a finite pull-off force, and as 
long as the contact area is small compared to the nomi¬ 
nal contact area. It is clear that if the nominal contact 
area in case (b) is very small, the latter assumption will 
no longer hold. But if all the assumptions made above 
holds, introducing a periodic “roughness” on the rubber 
surface will have no influence on the friction coefficient 
and the contact area. We will show in the next section 
that random roughness on the rubber surface does influ¬ 
ence the friction, but in most cases only to a very small 
extent. 

3 Theory: special case 

In the following we will extend the contact mechanics 
and rubber friction theory developed in Ref. @ to the 
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case where both solids have surface roughness, but one 
of those is rigid. In particular, in Sec. 3.1 we consider 
the contact area and in Sec. 3.2 the rubber friction. The 
most general case of arbitrary solids rheology is instead 
discussed in Sec. 5 and in Appendix A, where we present 
analytical results for the contact area, the rubber friction, 
and average interfacial separation. 


3.1 Contact area 


In this section we study the contact mechanics for the 
most important case where solid 1 is viscoelastic (say 
rubber) and solid 2 rigid, and both surfaces have surface 
roughness which are uncorrelated so that (/ii(x)/i 2 (x)) = 
(/ii(x))(/i 2 (x)) where (..) stands for ensemble averaging. 
In this case, in order for solid 2 to make perfect contact 
with solid 1, one needs to deform the surface of solid 1 
so that the normal surface displacement Uz (x, t) equals 
the gap function: 

■Uz(x,t) = /l 2 (x,t) - /li(x). (1) 

Let us define the Fourier transform 

Uz(q,a;) = J cPxdt Uz(x, 

If the upper solid moves at a constant velocity v = 
(vx,Vy) parallel to the bottom solid, we have /i 2 (x, t) = 
/i 2 (x — vt) and 

h 2 (q, w) = /i 2 (q)(5(a; - q • v) 

and ?? takes the form 

Uz(q,a;) =/i 2 (q)(y(a; - q • v) -/ii(q). (2) 

If a normal stress or pressure Uz (x, t) acts on the surface 
of solid 1, from the theory of viscoelasticity @ 

Uz(q,uj) = M(q,uj)az(q,uj). ( 3 ) 

In the Persson contact mechanics theory the mean-square 
fluctuation of the (normal) interfacial stress for complete 
contact (see Appendix A in Ref. i) is needed to quantify 
the diffusion function: 

(cr^(x,t)} = J (Pqduj |M(q,w)|"2 

x(uz(q,a;)uz(-q,-w)), (4) 

where Aq is the surface area and to the time over which 
we perform averaging. Using ?? we get 

{uz{q,uj)uzi-q,-uj)) = 

(/ii(q)tii(-q))[(I(a;)]^ -f (/i2(q)/i2(-q))[(5(w - q - v)]^ 
where we have used that /ii(q) and /i 2 (q) are (assumed 


to be) uncorrelated. Next, using that 

{hiiq)hi{-q)) = j^Ci{q) (5) 

(27r)^ 

(/i2(q)/i2(-q)) = 7^C2iq) 

{2Try 

where Ci (q) and C 2 (q) are the power spectra of the sur¬ 
faces of solids 1 and 2, respectively, and that 

[<5(a;)]2 = <5(a.);^ 


we get 


Aoto 


(uz(q,w)Mz(-q,-w)) = 
[C'i(q)5(a;) -I- C'2(q)(5(a; - q • v)]. 


( 6 ) 


(27r)3 

Substituting this result in ?? gives 

((T^(x,t)) = J d^qdu) |M(q,w)|“^ 

X [C'i(q)(5(w) -b C 2 {q)S{uj - q • v)] 

= J d^q [|A^(q,q-v)r^C2(q)-b |M(q,0)r2C'i(q)] . 


For a homogeneous viscoelastic solid 


M {q,uj) = -- 


1 E{uj) 


2 1-1/2 

Substituting this in ?? gives 

(CT^(x,t)) = i y q^ 
E{q ■ v) ' 


q- 


(7) 


( 8 ) 


I — :/2(q • v) 


C2{q) + 


m) 


I - :/2(0) 


c^i(q) • (9) 


If the rubber surface is smooth (A (*4) “ 9, ?? reduces 
to the expression derived in Ref. Q ■ For an elastic solid 
E(iij) = i?(0) and v{io) = i/(0) are frequency independent 
so in that case 


(cr^(x,<)) = ^Jd^q q^ 


E 


I - 1/2 


[C 2 (q) -bC'i(q)]. 


This result is the expected one because for a rigid solid in 
contact with an elastic solid the surface-mapping effective 
modulus E/2{l — iP‘) is just that of the elastic solid, while 
the effective power spectrum is that for the combined 
roughness given by 

C(q) = ■^^^([^i(q) + ^ 2 (q)] [hi{-q) + h 2 {-q)]), 
i.e. C'(q) = Cl (q) + C 2 (q), where we again have assumed 
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uncorrelated roughness profiles. Thus, for elastic solids 
the results above could be obtained directly using the 
surface-mapping. 

In the Persson contact mechanics theory the relative 
contact area is given by 

where erf is the error function, with the integral repre¬ 
sentation 


2 /■“ 

erf (a;) = —j— / dx exp(—a;^), 

VTr Jo 



and where 


x[|M(q,q-v)| 2 C' 2 (q)-h |M(q, 0)| ^C'i(q)] , (11) 


where the integral is over all the roughness wavevectors. 
The nominal contact pressure tro = Tq/Aq, where Fq is 
the normal force. For a homogeneous viscoelastic solid 
this expression takes the form 


G = 


8crg 




|F,(q-v)rC2(q) + |F,(0)rC'i(q) 


Figure 4. The surface roughness power spectra of an asphalt 
road surface and a rubber tread block surface after use. 

given by Eq. ?? in Ref. 

AE = j d^xdt 

= {2Trf J d^qduj {iu})M~^{q,uj){uz{-q,-Oj)uz{q,u;)). 

Next using ?? this equation takes the form 


where where (w) = E (w) / \l — v'^ (w)]. Note that G 
depends on the range of surface roughness included in the 
analysis. Thus, if go correspond to the smallest surface 
roughness wavevector, and if we include all the rough¬ 
ness qo < q < qi then the function G{qi) will mainly 
depend on the cut-off wavevector qi. In general, as gi 
increases, G{qi) increases and the contact area A{qi) de¬ 
creases. For the rubber friction to be discussed next we 
need the relative contact area P{q) = A{q)/AQ for all q 
between go and gi. 

Using that for a: << 1 we have erf(x) « 2x1^7^ we get 
for large G that AjA q k, Note that G —>■ oo 

as (To 0 so for small nominal contact pressures 
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/ l-£^r(q-v)| 

I \Erm 


2 

C^2(q) + Gi(q) 



AE = Aoto J d^qduj (*w)M“^(q, w) 

X [Gi(q)^(a;) -b G2(q)(5(a; - q • v)] 

= Aoto y d^g (iq-v)M“^(q, q • v)G 2 (q). (13) 

As expected, for complete contact the dissipated energy 
does not depend on the surface roughness on the rubber, 
but only on the roughness on the hard counter-surface. 
As a result the rubber friction coefficient within the Pers¬ 
son approach is given by the same expression as derived 
in Refs. Sin, except now enters the contact area func¬ 
tion P{q) given W ?? and ??, rather than the expression 
derived in Ref. Si a perfectly smooth rubber surface: 

d = ^ / dq q^C 2 {q)P{q)S{q) 

X f d(j) C0S(j) Im — f (14) 

Jo o-o(l - ^ ) 


It is clear that the surface roughness on the solid 1 will 
reduce the contact area. 


where P{q) = A(g)Mn is given by ?? and where the 
correction factor 


3.2 Rubber friction 

We now study the contribution to rubber friction from 
the viscoelastic deformations of the rubber surface in¬ 
duced by the surface ro^hness of the hard counter sur¬ 
face. Following Ref. Si ^^st focus on the case of 
complete contact. In this case the dissipated energy is 


F(g) = 7+(l-7)p2(g) 

with 7 Ri 1/2. Note that P(g) [and S'(g)] depend on the 
applied (or nominal) pressure (Tq and sometimes, when 
necessary, we write P{q) = P(g, (Tq). 

4 Numerical results 












V 



-12 -10 -8 -6 -4 -2 0 2 4 

logio V (m/s) 


Figure 5. The calculated contact area (a) and friction coef¬ 
ficient (b) for a rubber tread block sliding on a road surface 
with the surface roughness power spectrum shown in Fig. ?? 
(red curve). The red curve is assuming the rubber surface 
is completely flat (no surface roughness) and the blue curve 
assuming the same surface roughness on the rubber surface 
as on the road surface. If one instead use the actual surface 
roughness on the tread block (green line in Fig. ??) the calcu¬ 
lated friction coefficient and contact area becomes so similar 
to the case of the flat rubber surface that it overlap the red 
curves. For the nominal contact pressure (tq = 0.15 MPa and 
the temperature T = 20°C. Calculations are without includ¬ 
ing the effect of frictional heating. 


In this section we will illustrate the theory presented 
above with an example. Fig. ?? shows the surface rough¬ 
ness power spectra of an asphalt road surface and a rub¬ 
ber tread block surface after use. Note that, as expected, 
the road surface exhibits much larger roughness then the 
tread block. 

By using the theory above one finds that the friction 
(and the contact area) between the asphalt road surface 
and the rubber surface with the power spectra given in 
Fig. ?? are nearly identical to the case when the rubber 
surface is perfectly flat. For this reason, we instead as¬ 
sume that the rubber surface is as rough as the asphalt 
road surface i.e. it has a power spectrum given by the 
red curve in Fig. ??. For this case we show in Fig. ?? 
(blue curves) the calculated contact area (a) and friction 
coefficient (b). The red curves in the same figure are the 


results when the rubber surface is completely flat (no sur¬ 
face roughness). It is remarkable how small influence the 
surface roughness on the rubber block would have on the 
contact area and the friction even in this extreme case. 
Note in particular that for high sliding speeds the rough¬ 
ness on the rubber surface has negligible influence on the 
calculated friction. This is easy to understand from ?? 
and ??. Indeed, the contact area depends on the sum 



-Er(q-v) 
CTo 


2 

C'2(q) + 



. (15) 


For small sliding speeds, w = q - v is very small in a large 
part of the q-plane, and the w = q • v and a; = 0 terms 
in ?? are of nearly the same magnitude in a large part 
of the q-integration domain. However, for high sliding 
speeds w = q-v will be very large in a large part of the q- 
plane, and since E(ui) increases strongly as the frequency 
increases from the rubbery region (small w) to the glassy 
region (large w), it is clear that for large sliding speeds 
the w = q • V term in ?? will be much larger than the 
w = 0 term in a large region of the q-plane. This explain 
why the two curves in Fig. ?? (a) and (b) approach each 
other as the velocity v increases. 

5 Theory: general case 

We now consider the general case where both solids 
have surface roughness and arbitrary viscoelastic prop¬ 
erties. Using the Persson contact mechanics theory, in 
Appendix A we derive expressions for the contact area, 
the viscoelastic contribution to the sliding friction, and 
the average interfacial separation. In this section we con¬ 
sider 4 limiting cases (see Fig. ??) of the full theory. 

We first summarize the basic results obtained in Ap¬ 
pendix A. We assume that solid 1 moves with the velocity 
vi and solid 2 with the velocity V 2 . All physical proper¬ 
ties depends only on the velocity difference v = V 2 — Vi 
so we can always assume V 2 = v and vi = 0. The contact 
area is given by ?? with: 



Cj (q) 

\M (q,q-Vj) |2’ 


(16) 


where 

M (q,a;) = Mi (q, w - q • vi)-kM 2 (q,a; - q • V 2 ). (17) 

The viscoelastic contribution to the friction coefficient is 
given by: 


= 


<fqS{q)P{q)Y, 


t=i 


M (q, q ■ Vj) Cj (q) 
(Jo\M (q,q- Vj) |2 ’ 


(18) 


where 

M (q,a;) = Mi (q, w - q • Vi) -|- M 2 (q, w - q • V 2 ) 
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depends on [see ??]: 

A M (q,q-Vj-)Q (q) 

|M(q,q-Vj)|^ 

^ i (q- Vj - q - Vi)Cj (q) ^ iq ■ VC2 (q) 

Ml (-q, -q • (v^- - Vi)) Mi (-q, -q • v)' 

( 20 ) 


Note that ?? correspond to ?? as expected. 

Finally, the average interfacial separation depends on 
[see ??]: 

^ C,(q) Ci(q) C2(q) 

^M(q,q-Vj) Mi(q, 0 ) Mi(q,q-v)’ 


Figure 6. Four limiting cases discussed in the text. 


and 


Mj (q, w) = iojMj (q, w) . 

The average interfacial separation is given by: 

u{ao)=J d'^q'^Cj {q)Re M(q,q-Vj)”^ 

1 d[S{q,p)P{q,p)] 


dp 


' (To 


2p 


dp 


(19) 


Case I: A rigid solid 2 sliding on a viscoelastic solid 
1. Both solids have random surface roughness. 


As discussed in Sec. 4, we note that, at least for high 
sliding speed, the power spectrum of the compliant solid 
usually plays a negligible role in the generation of the 
contact area. Similar, the average interface separation u 
is determined mainly by the power spectrum of the rigid 
solid, but note the difference in the potency of the Mi- 
term (two for the contact area and one for the average 
separation). 

Case II: A rigid solid 2 with smooth surface sliding 
on a viscoelastic solid 1 with rough surface. 

With M2 = 0 and /12 = 0 , we have 

M (q,a;) = Mi (q,w - q - Vi) 

M (q, w) = i (w — q • Vi) M (q, w). 

The contact area depends on: 

Cj (q) _ Cl (q) 

^ |M(q,q-Vj)|^ |Mi(q, 0 )|^ 

The viscoelastic contribution to the friction depends 
on: 


This limiting case is the same as the theory developed 
in Sec. 3. With M 2 = 0, and considering that v = 
V 2 — vi, we have 


A M (q, q ■ Vj) Cj (q) 

j=i l^(q,q-vj)|^ 


M (q,w) = Ml (q,w - q - Vi) 

M (q, w) = i (w — q • Vi) M (q, w). 

The contact area depends on [see ??]: 

Cj (q) ^ Cl (q) ^ C2 (q) 

^.^1 |M(q,q-Vj)|^ |Mi(q, 0 )|^ |Mi(q,q-v)|^ 

Note that this equation correspond to ?? as expected. 
Similarly the viscoelastic contribution to the friction 


Hence, as expected, no hysteresis friction is generated 
during the sliding interaction. This is in accordance with 
the qualitative discussion presented in Sec. 2. 

The average interfacial separation depends on: 

(24) 

^M(q,q.v,) Mi(q,0)- ^ ’ 

Note that both the contact area and the average inter¬ 
facial separation depends on the viscoelastic modulus at 
vanishing frequency (i.e., in the rubbery regime). This is 
of course expected as there are no time-dependent defor¬ 
mations of the rubber surface. 



































i (q ■ v) Cl (q) 
M2(-q,q-v)' 
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Case III: An elastic solid 2 with a smooth surface 
sliding on a viscoelastic solid 1 with rough surface. 

With M 2 (q, w) = M 2 (q) and ^2 = 0, we have 

M (q, uj) = Ml (q, w - q • vi) + M2 (q, a; - q • V2) 

M (q, w) = Ml (q, w - q • vi) + M 2 (q) . 


on: 

^ M (q,q-Vj )Cj- (q) ^ 

Assuming v = V2 — vi is aligned along the x-axis, using 
?? in ?? we obtain 


The contact area depends on: 

^ Cj (q) ^_ Cl (q) _^ 

^ |M(q,q-Vj)l^ |Mi (q,0)+M2 (q)l^ 

The viscoelastic contribution to the friction depends 
on: 

^ M (q, q ■ Vj) Cj (q) ^ -i (q ■ v) M2 (q) Ci (q) 

,=i (q,q-Vj)|^ |Mi (q, 0 ) + M2 (q)l^ 

(26) 

The average interfacial separation depends on: 


which for a viscoelastic half space reads (removing sub¬ 
scripts for simplicity) 

h- ( cto ) = ^ y ^^9 qxqS (q) P (q) Ci (q) Im [E^ (q^v)] 

which corresponds to ??. 

The average interfacial separation depends on: 


Cj (q) ^ Cl (q) 

^M(q,q-Vj) M2(q,-q-v)' 


(30) 


Cj (q) _ Cl (q) _ 

^M(q, q-Vj) Ml (q, 0 ) -f M2 (q) ■ 

We note first in ?? and ?? that the equivalent interfacial 
compliance, related to the contact area as well as to the 
average interfacial separation, is determined by a simple 
summation rule of each solid compliance, evaluated un¬ 
der static interaction. Moreover, even if ?? is non-zero, 
(q • v) M 2 (q) is an odd function of q, resulting into a 
zero micro rolling friction [see ??]. 

Case IV: A viscoelastic solid 2 with a smooth surface 
sliding on a rigid solid 1 with rough surface. 

With Ml = 0 and ^2 = 0 we have 

M (q, w) = M2 (q, w - q • V2) 

M (q, w) = i (w - q • V2) M (q, w) . 

The contact area depends on: 

^ Cj (q) _ Cl (q) 

f^i\M{(l,q-Vj)f |M2(q,-q-v)|^ 

The viscoelastic contribution to the friction depends 


6 Summary and conclusion 

We have studied the contact mechanics and sliding fric¬ 
tion between a viscoelastic solid and a rigid solid where 
both objects have random uncorrelated roughness. We 
have shown that for a rubber-like material, where the 
Young’s modulus may increase by a factor ^ 1000 when 
going from the rubbery region to the glassy region, rough¬ 
ness on the rubber surface will in most cases have a rather 
small influence on the sliding (velocity v) contact. How¬ 
ever, for static contact, obtained in the limit u —>■ 0, the 
roughness on both surfaces is equally important. As a 
numerical application we considered the contact between 
a rubber tread block and a road surface, and showed that 
during sliding the surface roughness on the rubber sur¬ 
face has in most cases only a very small influence on the 
contact area and the viscoelastic sliding friction coeffi¬ 
cient. 

We have also studied the general case where both solids 
have random roughness and arbitrary (linear) rheology 
properties. For this case, within the Persson contact me¬ 
chanics approach, we have derived expressions for the 
contact area, the viscoelastic contribution to the friction 
force, and for the average surface separation. We have 
studied 4 different limiting cases of the theory. 

The theory presented in this paper can be extended 
to correlated surface roughness, and to include frictional 
heating and adhesion. 
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Appendix A: General contact mechanics theory 


Here we derive the general contact mechanics theory for steady sliding rough solids (to which a moving reference 
is associated, where i = 1 for the lower solid, see Fig. ??), in adhesive- and friction-less interaction. We define 


Mz (q,^,w) 



ix,z,t) 




(Al) 


where Uz is the generic surface out-of-plane displacement, and 


cr(q,z,a;) 



z,t) 


(A2) 


where cr (x, 2 ; = 0, t) = cr (x, t) is the contact stress, both observed in the fixed reference Hence, given that 
[resulting into cr|^ (q, w) = (q, w — q • v^), and similarly for Uz], we get 

cr (q, w) = M~^ (q, w - q • v^-) Uzj (q, w), (A3) 


whereas in term of the displacement time derivatives [note: [uzj (q, w)] 

resulting into z (u; — q • Vj) Uzj (q, to) = Uzj (q, u;)] 




cr (q, w) = MO - q - Vj) Wj (q, w), 


where simply 


Mj (q, w) = iu>Mj (q, w). 


Hence, by equating the stress in ?? for both top and bottom solid we get 

Uz2 (q, w) = M2 (q, w - q • V2) Mf ^ (q, w - q • vi) Mzi (q, w) 


so that 


Uz (q, w) = Uzi (q, w) + Uz2 (q, w) = [l + M2 (q, w - q • V2) M;^ ^ (q, w - q • vi)] Mzi (q, w) . 
We can determine the equivalent interface compliance M (q, w) from 

cr (q, w) = M"i (q, w) Uz (q, w), 


(A4) 

(A5) 

(A6) 

(A7) 


where 


M (q, w) = Ml (q, w - q • vi) -I- M2 (q, w - q • V2). (A 8 ) 

Moreover since iizi (q, w) /Mi (q, w — q • v^) = Uz (q, w) /M (q, w) we can determine the equivalent interface compliance 
derivative M (q, w) from 

cr (q, w) = M"i (q, w) iiz (q, w), 

where 


M (q, w) = Ml (q, w - q • vi) -I- M2 (q, w - q • V2) . 
Observe finally that from ?? and ?? it follows 


M ^(q,a;)M ^(-q,-w) 
M-i (q,w)M-^ (-q,-w) 


|m (q,w)| 
|M (q,w)|"^ 


i.e. [M ^(—q,— w) = [M ^(q,a;) * [and similalry for M ^(q,a;)]. 


(A9) 



IX 


We calculate first the stress auto-correlation for full contact. In such a case 

2 


u. 


(q, w) = ^ hj (q) 5 (w - q • Vj) 
i=i 


i.e. 


CT (q, w) = M ^ (q, w) ^ hj (q) 5 (w - q • v^). 

Hence, considering that Cj (q) 5 (q + q) = {hj (q) hj (q)) and assuming the surface roughness are uncorrelated 


(ct^(x, t))= / (f'qdw I d^qdoj (cr (q, w) cr (q, w)) ‘^*)g*(q-x <^t) 


= /ae- 


Q (q) 


j=i\M iq,q-Vj)\ 

With similar considerations we can calculate the correlation function (a (x, t) (x, t)) 


(AlO) 


(cT (x, t) Uz (x, t)) = J d^qduj J d^qduj (a (q, cc) Uz (q, a;)) e 
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z(q-x—u;t) (q-x—u;t) 


= J d\S{q)P{q)Y^ 


M (q, q • Vj) Cj (q) 


|A^(q,q-Vj)| 

Furthermore, the correlation function {a (x, t) (x, t)) reads 

(ct (x, <) Mz (x, t)) = / d^qdw [ d^qdu) (ct (q, w) Wz (q, w)) 


(All) 


= J d'^q S (q) P {q) E Q (q) Re [M (q, q • v^)]. 


(A12) 


In ?? and ?? we have added the S {q) P (q) term to take into account the partial contact in the evaluation of 
(ct (x, t) ftz (x, t)) and (cr (x, t) Uz (x, t)), where S (q) has been introduced before. Mi{q,ui) = 2 (w) S (qd)] 

for a generic viscoelastic solid of finite thickness, where the real function S{qd) is a correction factor related to the 
adopted boundary conditions [e.g. for a slab of thickness d on a rigid substrate we have 


S{qd) = 


(3 — Ah') cosh {2qd) + 2 (qd)^ — Ah' {3 — 2i/) + 5 


(3 — Av) sinh {2qd) — 2qd ’ 

where S (qd) —>■ 1 for d—>■ oo]. 

Contact area 

The contact area A/Aq = P (qi) can be calculated from ??, with G (q) = {a^ (x, t)),/ (2crQ) given by ?? 


— 2 

G[h) = i / AE 


C, (q) 


^ |Ar(q,q-Vj)| 


(A13) 


Stored and dissipated interfacial energy 

In the process of determining the average interfacial separation u as well as the micro-rolling friction, we assume the 
normal approach to occur quasi-statically with respect to the sliding kinematics. Within this assumption (i.e. to tij 
where to (x u~^ and ti (x v~^ are, respectively, the time scale of the normal and sliding motion), the time-derivative 
of the displacement field iiz (x, t) can be expanded into two uncorrelated sources wq (x, t) + wi (x, t) (ivi is random in 



X 


time), where uiq (x, t) cx tg ^ and wi (x, t) oc Here we calculate the average external work over a ti time span as 


W^= dtW= dt [—aoAou + toAqv] . 

Jti Jti 

Since ti/to ^ 1 the term agAgii can be assumed constant and given that dt tq = (r), we get VH, = —ticrgHott + 
tiAov (r). However u oc tg ^ and v oc t^^, resulting into 

ft! hAov (t) . 

Now the average external work over a tg time span reads 


W = f dt W = f dt [—(TgHgu] + — W^. 
Jto Jto tl 


(A14) 


However the energy conservation requires that W = Aq (a (x, t) w (x, t)} tg = Ag (tJo (x, t) wg (x, t)} tg + 
Ag {ai (x, t) ibi (x,t)} tg. However since (dg (x,t) wg (x,t)) oc tg^ and (cti {x,t)ibi (x, t)) oc t^^ we can separate the 
different energy contribution in ?? as 


(ag (x, t) Wg (x, t)} tg= dt [-CTom] 
dto 


(A15) 


Ag (cTi (X, t) Wi (X, t)) ti=W^. 


Manipulating the LHS of ?? we have 


(cto (x, t) Wg (x, t)) tg = (ag (x, t) Awg (x, t)) = - A (ag (x, t) Wg (x, t)) , 


resulting into 


(A16) 


-A (ag (x, t) Wg (x, t)) = -agAu, 

whereas ?? leads to 

(cti (x, t) wi (x, t)) = (v 2 - vi) • (r). 

Hence, by using ?? in ?? we get the sliding-induced micro-rolling shear stress 


r ^ 

(V 2 - vi) • r (p) = J d?q S (q) P (q) ^ 


^ (q, q • Vj) Cj (q) 


|2 ’ 


j=i |A4^(q,q-vj)r 

whereas the average interfacial separation can be calculated by substituting ?? into ??, 

/» 2 /.on 


'.(p)= / d^q Cj (q) Re (q, q • v^ ) 
i=i 


-1 




1 d[S{q,ao)P{q,ag)] 


dai 


0 


(A17) 


(A18) 


(A19) 


(A20) 


For an elastic solid in contact with a hard randomly rough surface, this equation reduces to the equation derived in 
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